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Abstract 

We present some model calculations of persistent currents in dis- 
ordered one- and two-dimensional mesoscopic systems. We use the 
tight-binding model and calculate numerically the currents in small 
systems for several values of disorder. Next we fit appropriate ana- 
lytical formulae, and using them we find self-sustaining currents and 
critical fields in larger, more realistic systems with different shapes of 
the Fermi surfaces (FS). 
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1 Introduction 



With advances in technology, fabrication of submicron devices has become 
possible. Such systems exhibit quantum coherence which is the subject of 
extensive experimental and theoretical studies. 

Recently, in a series of papers [1], model considerations of self-sustaining 
persistent currents in a set of clean mesoscopic metallic rings have been 
presented. In this paper we investigate whether such currents survive in 
the presence of scattering driven by impurities and temperature. We will 
consider a set of concentric mesoscopic rings stacked along a certain axis and 
a set of concentric two-dimensional (2-D) cylinders. Such samples can be 
obtained by lithographic methods from metals or semiconductors 0. 

2 Quasi One-dimensional Mesoscopic Rings 

Let us consider a set of quasi one- dimensional (1-D) metallic mesoscopic rings 
stacked along y axis threaded by the external magnetic flux <fr e . We assume 
that each ring posseses the same number of conducting electrons N e (even or 
odd). The situation when some of the rings carry an odd number and other 
carry an even number of electrons has been considered in [|la[] . The long range 
magnetostatic (current-current) interaction between electrons from different 
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rings is taken into account in the mean field approximation (MFA). 

The gauge-invariant tight-binding Hamiltonian for a single ring is of the 
form: 

N 

% =]T [{2t + V n )c+c n - te ie "^ct +1 c n - te- ie -" +1 c+c n+1 ], (1) 

71=1 

where t = ft 2 / (2m e a 2 ) is the hopping matrix element, a is the lattice constant; 
c+,c n are the creation and annihilation operators; N is the number of sites 
in each ring (channel); 6 n , n +i comes from the magnetic flux <p : 

^ = (Ja-«=%± = 0, (2) 

where A is the vector potential, and 0o — h/e. 
The flux (f) is composed of two parts, 

= e + 7 , 0j = £/(0), (3) 

i.e., each electron moves in the external magnetic flux <p e and in the flux com- 
ing from the total current / in the system; £ is the selfinductance coefficient. 

The disorder is given by a random choice of the on-site potentials V n from 
a rectangular distribution of width W — nt, k > 0, and strength — W/2 to 
W/2. 

In the case of a clean system (V n = 0) we can calculate the current I in 
the system, diagonalizing the Hamiltonian ([!]) directly, by using the Fourier 
transform. We get 

7(0) = M t h{<f>), 
3 



= - Yl if = I — 2 /« sin 



2tt 



where £ n = 2t 
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(4) 

f a is the Fermi-Dirac distribution func- 
tion, a is the orbital quantum number for an electron going around the ring 
(a = 0, ±1, ±2, ...); M t = MP, M and P are the numbers of rings along and 
perpendicular to y axis respectively (P <C N). 

In this paper we concentrate on a disordered system (V n ^ 0). The 

Hamiltonian (|l|) can be written in the form of a matrix (NxN) : 

/ 2t + Vi -te~ iei ' 2 -te i0N ^ 

-te idl - 2 2t + V 2 -te~ i62 - 3 

— -te ie ^ 3 2t + V 3 -te'^ 4 

Hi = 

_ te ie N . 2 ,N-i 2t + V N -i -te~ i0lf - 1 ' lf 

V -te-^- 1 _ te ie N . UN 2t + V N J 

(5) 

where the base vectors are chosen as follows: 

c+ |0) = |100...0) = |1) 
4 10) = |010...0) = |2) 



c+ |0> = |000...1) = I AT) 
|0) = |000...0) - the vacuum state. 

The above matrix we have diagonalized numerically for different N (<100) 
and different W/t using the "Monte Carlo" method |4[]. We also calculated 
the current I\{4>, W/t, T) by numerical differentiation; it is presented in Fig. 
1 for iV = 80 by symbols ||. We see that the impurity scattering leads to 
the amplitude reduction and the shift of the maximum of the Ii(<f>) charac- 
teristics. We note M that elastic scattering has a very similar effect on the 



Ii(<f>) as temperature. Whereas a non-zero temperature leads to a redistribu- 
tion of electrons among the energy levels, collisions lead to a change in the 
levels themselves. If the average collision time is r, then the uncertainty in 
the electron energy is of the order of h/r, inverse proportional to the elastic 
scattering mean free path l e ~ (t/W) 2 || (the localization length £ in one 
dimension). If r is not too short, it plays a similar role as the temperature 
1- 

Because of the limited capacity of computers and long time needed for 
computing, it is very difficult to calculate in this way the current h((f>, W/t, T) 
for N > 100. However we want to investigate the possibility of the existance 
of self-sustaining current in the system with impurities. The phenomenon of 
such current is a collective effect ||1 b|| and can be obtained only if the number 



of interacting electrons is relatively large, i.e. for large (but still mesoscopic) 
systems. The selfconsistency formula for the current at <p e = 0, obtained by 
the use of Eq. (|||), is of the form: 

I^,W/t,T) = ^, (6) 

where 

C=^(VP^W-R), (7) 

/ = Mb is the length of the cylinder made of a set of mesoscopic rings, b is 
the distance between rings; fio is the magnetic permeability of free space. 
To get a stable, nontrivial solution of Eq. (0) we need to consider a set 



of mesoscopic rings with N ^> 100. Even though exact numerical solutions 

of the current in the presence of impurities were performed for very small 

samples, they show generic features for all mesoscopic rings [[|. Thus we 

have fitted an analytical formula to the numerical results of Fig. 1 and we 

assumed that it can be used for the qualitative analysis of Eq. (P) for larger 

rings. Ji is periodic in 0/0o? with period 1, and can be expressed as a Fourier 

sum The best fitting, denoted in Fig. 1 by solid line, we have obtained 

for the formula: 

t/i ^ AIq ( L T\ exp \-q (± + 

^ vr \27 T*J i_ exp [- g (| + 2^)J (g) 

x cos(qk F L) sin f 27rg— ] , 
V V 

where Iq = eh sm(kpa) / (Nm e a 2 ) , kp = N e ir/L, L = Na is the circumference 
of a ring; T* = A/(27t 2 /cb) is a characteristic temperature that separates the 
high- and low-temperature regimes, A = 4nt sin(kpa)/N is the level spacing 
at zero flux at the Fermi surface (FS), ks is the Boltzmann constant; I/7 is 
a disorder parameter, 7 = 10a(t/W) 3 ^ 2 In N. 

To account for different values of 7 (W/t) in a set of M t mesoscopic 
rings we have performed a quenched average of the current I ((f), 7, T), where 
disorder was given by a Gaussian distribution with standard deviation A7 = 
±10%7, 7 is a mean value of 7. We have obtained a very small decrease 
in the amplitude of [7(0,7, T)], where [...] denotes the quenched average, 
comparing to the case with the same disorder in each ring. 
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In the clean system at T = K we would have ||: 



oo 



21, 



(po 



J(0) = M t £ 



o 



cos(qkpL) sin 27rg 



(9) 



Now we are in position to look for the self-sustaining currents. They can 
be obtained from the self-consistent equation for the total current (^), where 
the current / = M t l\ is taken with the quenched average, I\ is given by the 
formula (Q). 

The case of a clean system, in the free electron model, has been analysed 
in ||lb|| . Here we extend earlier studies to include the dependences on disorder 
and temperature in the tight-binding model with the half-filled band. The 
results depend on the number of electrons in a single ring. 

In the case of an even number of electrons in each ring the current is 
paramagnetic and the stable solutions of Eq. (|GJ) being the intersections of 
the curves marked by circle correspond to a spontaneous current I s . They 
are presented in Fig. 2 for different parameters 7. In Fig. 3 we also present 
the temperature dependence of persistent current and we see that both tem- 
perature and disorder act almost in the same way. 

In the case of an odd number of electrons in each ring the current is 
diamagnetic for small and the stable solutions of Eq. (Q) presented in Fig. 
4 correspond to flux trapping, the phenomenon known in superconductors. 

The temperature at which the transition to the state with the self-sustaining 
current occurs is denoted by T c . We see that T c decreases with increasing 
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the impurity content and that flux trapping can be obtained only in weakly 
disordered systems. 

3 Two-dimensional Mesoscopic Cylinders 

Let us consider now a set of coaxial closely packed two-dimensional (2-D) 
cylinders [fLl]] . The appropriate tight-binding Hamiltonian for a single cylin- 
der is of the form: 

_ M N 

H = 2^ [(2i + Vnm) c nm C nm ~ te n,n+1 C n _|_ lm C nm + 

m=lra=l (^10) 

t€ n ' n+1 C nm C n j t .\ m tf 1 C nm j r iC nm tfeC nm C nm ,_|_i] , 

where tb = h 2 / (2m e b 2 ) is the hopping element in the y direction, b is the 
lattice constant in the y direction. 

In this case to get the current J in a disordered system, we have performed 
similar procedures to the previous 1-D case. Finally we have obtained: 



J(0, 7 ,T) = P£E 

m=lq=l 



Vl x > 4J (m) / L T 
n \ 27 + T* 



exp 


[-*(*+£) 






1 — exp 







X 



x cos[qkF x {m)L] sin 27rg— , 
V 9o J 

(11) 

where I/7 is a disorder parameter, 7 ~ (t/VT) 3 / 2 In N; 

Io(m) = ehsm[kF x (m)a]/(Nm e a 2 ), kp^irn) is the radius of the FS corre- 
sponding to the channel m which depends on the shape of the FS. 

After performing a quenched average of 7(0, 7, T) over 7 iW/t) with the 
Gaussian distribution we can look for self-sustaining currents. 



In the clean system at T = K we would have again ||: 

I{<P)=P EE ^^cos(g^(m)L)sinf2vrg^ . (12) 
m =i q =i *Q V 90 J 

The self-sustaining currents can be obtained from equations (^) and (|TT1) . 

It is known that for 2-D systems and for 3-D systems with 2-D conduction 

(an example of such structures are high T c superconductors in a normal 

state flQfl ) persistent currents depend crucially on the strength of the phase 



correlation between currents of different channels ||, |T0| and therefore we 
have analysed three different cases namely the systems with rectangular, 
triangular and half-circular FS. The k Fx (m) is then given by 

_ ±1 for \k Fy (m)\ < f, 
for other k F (m), 



k F (m) = ±— arccos 
a 



a 2 



b2 cos(k Fy (m)b) 



(14) 



and 



k Fx (m) = ±J^ + ^-k Fy (m) (15) 



respectively; k Fy {m) = mir/((M+l)b). Such shapes of the FS can be obtained 
in the tight-binding approximation depending on the crystal symmetry and 
on the band filling. 

The most favorable situation is for the rectangular FS (found in bcc 
crystals for nearly half filled band) where self-sustaining currents I s are the 
easiest to obtain and are the largest. In the second case (crystals with simple 
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cubic symmetry at half filling), it is possible to obtain self-sustaining currents 
only for weakly disordered systems (ballistic regime) and for the maximal 
interchannel phase correlation, where L/l = 2 + h/3, h is a positive integer; 
Fig. 5. In the third case (found at very low band filling), we have not obtained 
such currents, for physical values of parameters, because the interchannel 



correlations are too small | lQfl . The [I(4>)} characteristics for different shapes 
of the FS and for 7 = 50000 A are presented in Fig. 6. 

Self-sustaining currents are a hallmark of a phase coherence. Thus we 
can introduce, in analogy to superconductors, a notion of a critical field H c 
below which the system is in a coherent state. The system is in a coherent 
state if we get at <p e = the nonzero crossing of the currents J(0) and 4 
(Eq. (||)). Let us denote by max the value of <fi at which 7(0) reaches its 
maximum value J max = J(0 max ). 

If 

Z(</W) > , (16) 



we get a non-zero solution. Accepting the interpretation JT2[ that the critical 
field H c in the cylinder is determined by the maximum value of the current 
J max one finds from Eq. ([16]): 

H c = ^%|. (17) 

The magnitude of J max depends on temperature, disorder and on the geom- 
etry of the FS and can be extracted from Figs. 4, 5. 
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The critical field H c as a function of the disorder parameter I/7 and tem- 
perature T is presented in Figs. 7 and 8 respectively. We see that it decreases 
with the increasing disorder and likewise decreases with the increasing tem- 
perature. 

4 Conclusions 

In the presented paper we have addressed the question of whether self- 
sustaining currents in nonsuperconducting rings and cylinders survive in the 
presence of disorder and finite temperature. The answer to this question 
is positive. We have performed some model considerations which show that 
self-sustaining currents can be obtained in weakly disordered systems at tem- 
peratures of the order of 1 K in the systems with strong correlations among 
the channel currents, i.e. in the systems with the FS having flat regions. 

From the presented /(</>) characteristics we were able to determine the 
critical field H C) below which the system is in a phase coherent state. We 
have shown that H c decreases with increasing disorder and temperature re- 
spectively. 
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Figure Captions 

Fig. 1. Persistent currents lj in a mesoscopic ring as a function of 
magnetic flux 0/0o for different values of disorder. Numerical results for 
different W/t are presented by symbols and analytical results according to 
formula (§) for different 7 are plotted by solid lines. 

Fig. 2. Persistent currents as a function of magnetic flux 0/0o f° r 

different parameters 7, in a set of 1-D mesoscopic rings with even number of 
conducting electrons N e in each ring. Self-sustaining spontaneous currents 
I s are denoted by circles. 

Fig. 3. Persistent currents [I]/Iq as a function of magnetic flux 0/0o for 
different temperatures T, for the case presented in Fig. 2 with 7 = 20000 A. 

Fig. 4. Persistent currents [I]/Iq as a function of magnetic flux 0/0o for 
different parameters 7, in a set of 1-D mesoscopic rings with odd number of 
conducting electrons N e in each ring. Self-sustaining currents I s denoted by 
circles correspond to flux trapping. 

Fig. 5. Persistent currents [I]/Iq as a function of magnetic flux 0/0o for 
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different parameters 7, in a set of 2-D concentric mesoscopic cylinders with 
the triangular Fermi surface. Self-sustaining currents I s (trapped fluxes) are 
denoted by circles. 

Fig. 6. Persistent currents [/]/% as a function of magnetic flux 4>/4>o for 
different shapes of the 2-D Fermi surfaces for 7 = 50000 A. In the inserted 
figure the shapes of the FS are shown. 

Fig. 7. Critical fields H c as a function of the disorder parameter 7 _1 for 
the cases presented in figures 4, 5. 

Fig. 8. Critical fields H c as a function of the temperature T for the cases 
presented in figures 4, 5. 
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